Signatures of Fermi surface reconstruction in Raman spectra of underdoped cuprates 
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We have calculated the Raman Bi 9 and B2 9 spectra as a function of temperature, as well as doping, 
for the underdoped cuprates, using a model based on the resonating valence-bond spin-liquid. We 
discuss changes in intensity and peak position brought about by the presence of a pseudogap and 
the implied Fermi surface reconstruction, which are elements of this model. Signatures of Fermi 
surface reconstruction are evident as a sharp rise in the doping dependence of the antinodal to 
nodal peak ratio which occurs below the quantum critical point. The temperature dependence of 
the Bi 9 polarization can be used to determine if the superconducting gap is limited to the Fermi 
pocket, as seen in angle resolved photoemission spectroscopy, or extends beyond. We find that the 
slope of the linear low energy B2 9 spectrum maintains its usual d-wave form, but with an effective 
gap which reflects the gap amplitude projected on the Fermi pocket. Our calculations capture the 
main qualitative features revealed in the extensive data set available on the HgBa2Cu04+s (Hg-1201) 
cuprate. 

PACS numbers: 74.72.-h, 74.20.Mn, 74.25.Gz, 78.30.-j 



I. INTRODUCTION 

The underdoped region of the high T c cuprate phase 
diagram is known to contain a number of features which 
makes both the normal and superconducting state chal- 
lenging to understand using microscopic models. These 
features included: nanoscale inhomogeneities as seen in 
scanning tunneling microscopy (STM)^ the observation 
of a pseudogap energy scale above T,.^ as well as a 
transition to an antiferromagnetic Mott insulating state 
at low doping £ 

The discovery of a pseudogap feature in the normal 
state above the superconducting dome, for doping be- 
low some critical value (thought to be a quantum critical 
point), has recently resulted in considerable research ac- 
tivity. Understanding the effect of this pseudogap energy 
scale on the transition temperatures, T c , of the cuprates, 
as well as the origin of this energy scale itself, are impor- 
tant and interesting issues. Examination of the cuprate 
phase diagram shows a clear reduction in T c as the pseu- 
dogap energy scale increases. This could be an indication 
of the pseudogap arising from a competing phase, driv- 
ing down the superconducting gap, and hence the T c . 
Indeed, a great deal of experimental evidence appears 
to suggest that the superconducting gap and pseudogap 
occupy unique regions of phase space (nodal supercon- 
ductivity and antinodal pseudogap) iS 

Since its appearance in 2006, the model of Yang, Zhang 
and Rice (YRZ)- has shown some success in understand- 
ing Raman and optical properties'^ angle-resolved pho- 
toemission (ARPES)iS and specific heat data 1 - 1 . The 
YRZ model provides an ansatz for the coherent part of 
the many-body Green's function for an resonating va- 
lence bond (RVB) spin liquid system. In the absence of 
a pseudogap, the model maintains a large tight-binding 



Fermi surface that can undergo a transition to the su- 
perconducting state resulting in a superconducting gap 
which opens on this Fermi surface. The inclusion of a 
pseudogap in the YRZ model is such that it opens about 
the antiferromagnetic Brillouin zone (AFBZ) boundary. 
The net result is that a finite pseudogap in this model 
acts to deform the tight-binding Fermi surface to form 
Fermi hole pockets. If one assumes the onset of the pseu- 
dogap to be a zero temperature quantum critical point 
(QCP) at a critical doping x c , then we can define several 
regions around this critical doping. First, for dopings 
above x c , there is a superconductivity-dominated region 
with a large tight-binding Fermi surface. Second, for dop- 
ings well below x c , there is a pseudogap-dominated region 
with Fermi pockets which become smaller and shrink to- 
wards the nodal direction as the antiferromagnetic Mott 
insulator is approached. Finally, there is an intermedi- 
ate region wherein the features of a pseudogap onset, 
as well as a deformation of Fermi surface both present 
themselves. It is in this region that the energy scale of 
superconductivity and the pseudogap are comparable. 

There are other theoretical methods for including 
a pseudogap phase through, for example, preformed 
pairs' or as a competing d-density wave orderji 3 - As 
the YRZ theory has thus far resulted in good qualita- 
tive descriptions of a number of experimental properties, 
continued investigation of this model allows us to distin- 
guish quintessential features of the cuprates which must 
be included in any more complex microscopic models that 
might follow. Of central importance to this paper is the 
presence of two energy scales, nodal and antinodal, in the 
superconducting underdoped cuprates, as seen in Raman 
and angle-resolved photoemission spectroscopy (ARPES) 
experiments ! 14 ' 15 i 16 i 17 Within the model of YRZ, which 
includes two independent energy scales (superconducting 
and pseudogap), we calculate the minima in the angu- 
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lar band energies for comparison with ARPES and ap- 
plication to Raman spectra. Although the YRZ model 
has previously been applied to the calculation of Raman 
spectra at zero temperature in the work of Valenzuela et 
alf^, the work presented here extends upon those doping 
dependent calculations to understand better the specific 
behaviours which are present in models of the pseudogap 
phase. Here we will identify the impact on Raman Bi ff 
(antinodal) and Pj2 g (nodal) spectra of two important be- 
haviours: (1) the fc-space separation of superconducting 
and pseudogap energies and (2) the reconstruction of the 
Fermi surface. We will demonstrate how (1) and (2) play 
a role in the temperature and doping dependence of the 
Raman spectra. 

We have organized this work to fit the following struc- 
ture. Section II will describe the theoretical framework 
and parameters involved in the YRZ model used for all 
calculations shown. Section III will be a description of 
ARPES results in the context of the YRZ model with 
the motivation that this will be necessary for the under- 
standing of the Raman results. Section IV will contain 
our Raman calculations in the YRZ model, including fi- 
nite temperature, while Section V will contain a summary 
of our main conclusions. 



II. THEORY 

In the YRZ model£±£ both the superconducting gap, 
A sc , and the pseudogap, A pg , have a d-wave fc-space de- 
pendence described by: 



A° c (z) 



(cos k x a — cos k y a) , 



Apg = — — (cosfc^a — cos k y a), 



(1) 
(2) 



where a is the lattice constant. For a doping x, the YRZ 
model employs, for the coherent piece, a propagator, 



(3) 



which has been formulated from numerical studies of 
RVB-type models and will contain the essential physics 
of the superconducting state. Entering this Green's func- 
tion are the quantities: 

t _ to 



1 ± 



Ek 



(4) 



The energy dispersion = — 2t(cosk x a + cos k y a) — 
it' cos k x a cos k y a — 2t" (cos2k x a + cos2fc„a) — fi p is the 
third nearest-neighbor tight-binding energy dispersion, 
while £2 = — 2t(cos k x a + cos k y a) is the first nearest- 
neighbor term, which effectively shifts the placement of 
the pseudogap off the Fermi surface to an energy which 
coincides with the AFBZ boundary, defined by the half- 
filling point of the tight-binding energy ££. These en- 
ergy dispersions contain doping dependent coefficients: 
t(x) = g t {x)t Q + 3g s (x)Jx/8, t'[x) = g t [x)t' , andt"(x) = 



g t (x)t'o, where g t (x) = jfz and g s (x) 



JT+W' 



are the 



energy renormalizing Gutzwiller factors for the kinetic 
and spin terms, respectively. gt{x) also appears in Eq. © 
as a weighting factor for the coherent part of the Green's 
function which acts to statistically remove or project 
out doubly occupied states i 18 ' 19 The dispersion, uses 
/i p as a chemical potential determined by the Luttingcr 
sum rule.— Values of other parameters in the dispersion 
were taken from Rehl to be: t'/t = -0.3, t"/t = 0.2, 
J/t = 1/3, and x = 0.338 which are accepted values for 
the hole doped cuprates^I. The optimal superconducting 
gap Ag C was chosen to give an optimal T c around 95K 
for a ratio of 2A° c (x, T = 0)/k B T c = 6, using of t = 175 
meV. We also take Ag C (x, T), the gap amplitude at dop- 
ing, x, and temperature, T, to have a BCS temperature 
dependence. We will assume that the pseudogap is only a 
function of doping, and hold the value of A pg (a;) constant 
with temperature. This will allow us to attribute any 
temperature dependence in calculated quantities as be- 
ing due to the superconducting state. A° c (a;) and Ap g (a;) 
are described by the well known superconducting dome 
and pseudogap line, the latter vanishing at T = at a 
QCP in this model. These are given explicitly as 

A^ c (a;, T = Q) = 0.14t (l - 82.6(a; - 0.16) 2 ), (5) 
A° pg (x)=3t (0.2-x). (6) 

From the YRZ Green's function of Eq. ([3]) and stan- 
dard equations of superconductivity for the anomalous 
propogator, one can extract the regular and anomalous 
spectral functions, A(k,uj) and B{k 1 uj) respectively, and 
see that there are four energy branches, given by the en- 
ergies, 



±E*=±y/ + (7) 
These energy branches appear in the spectral functions 



A(k, w) = £ gtWZ[{u a )H{Lo - E?) + (v a ?8{Lo + £?)], 

(8) 

B(fc,w) = 9tWZu a v a [6(Lo + E«) - 5{u - E«% 

(9) 



a=± 
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where 
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E" 



»\iV 2 



(10) 

(11) 



The Raman response, = I m \Xv(Ti is given by 

poo 

x [A(fc, w)A(fc, w + n) + B(k, uj)B{k, cj + n)], (12) 

where rj is the choice of vertex Bi g or B2 5 . The ver- 
tex strength, (j^) 2 , can be determined straightforwardly 
from the energy dispersion, such that 



7 fe 
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,,,,, _ a 2 e(fc) 



dkx dky 

For the energy dispersion used, this results in 

Bi g _ 2ta 2 (cos k x a — cos k v a) 
+ 8t"a 2 (cos2k x a — cos2k y a), 
"/ k 23 = — At' a 2 sin k x a sin k y a. 



(13) 
(14) 



(15) 
(16) 



The Raman response can be simplified by performing 
the integration over to in Eq. (TT2"]) . This allows us to sep- 
arate the contribution to the Raman spectra from each 
energy branch so that 



X 



W - j E (rt) 2 9t{wz[f(E«)-f(E« + n)} 



h.a.— ± 



x [(u a ) 2 A(k, Ef + + u a v a B{k, Ef + ft)] 

+ W£\f(-E?)-f(-E? + n)] 

x [(v a ) 2 A(k, -Ef + fi) - u a v a B{k, -Ef + O)] 



(17) 



In the clean limit, Eq. (|17p can be written in the form 



x \u + u — v + v } 2 



6(Sl + B+ - E-) - S(uj + E- - E+)] 
+ [/(-£+)- f(E-)]W+W k 
x [v + u~ + v-u+} 2 S(n - E+ - E~) 
+ If(-Et) f(Er)](W+) 2 [2(u + v+) 2 }S(n - 2E+) 
+ [f(-E-) - f(E^)}(W k ) 2 [2(u-v-) 2 }S(n - 2E~) 

+ {HEt)(W+) 2 [{u+) 2 - [v+) 2 ] 2 

+h{E-){W k ) 2 [{u-) 2 - (v-) 2 ] 2 }n5(Q)) , (18) 




FIG. 1: (Color online) Plot of EJ band as a function of length, 
L, along a line traced from the (n,n) (see inset) point for a 
range of angles, 6, in the x = 0.14 doping case. The local 
minima in this curve are the energies of nearest approach, 
A na . The inset shows the fc-space locations of these near- 
est approach energies, marking out the Fermi pocket and its 
extension to the Brillouin zone boundary through the Fermi- 
Liquid-like 'wings', shown in red. 



where the function h(E) = — ■ In this form one 
can see the interband (first and second terms) , intraband 
(third and fourth terms) and Drude (last term which 
vanishes in the clean limit) contributions to the Raman 
spectra. However, with moderate impurities Eq. (| 18(1 be- 
comes inaccurate at low frequencies, due to not prop- 
erly capturing the Drude contribution. For this reason, 
throughout this work, all Raman calculations shown are 
the results of Eq. flT7|). 



III. ARPES 

In this section, we wish to understand how the angular 
energy profiles seen in ARPES data come about in the 
YRZ model. This will be essential to the Raman analysis 
and discussion given in Sections IV and V. 

The concept of a Fermi surface, defined here by = 0, 
can be a powerful aid in understanding the fc-space dis- 
tribution of electronic states and excitations. In the 
presence of energy gaps, more complicated energy dis- 
persions, E^, arise which may not trace a zero energy 
surface. In these cases, the examination of the minima 
of the energy dispersion produces a surface which rep- 
resents the lowest energy excitations in fc-space. In the 
YRZ model, examination of the E+ and E~ energies as 
a function of angle, serve to produce these lowest energy 
surfaces, or nearest approach energy, A na , surfaces. In 
Fig. [TJ the E~ energy is plotted as a function of dis- 
tance along a line traced from (7r,7r) for a range of an- 
gles. At 6 = 45°, the line is traced through the nodal 
direction, where both the pseudogap and superconduct- 
ing gap equal zero. This results in a true E~ = Fermi 
surface in the (7r,7r) direction, culminating at well known 
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FIG. 2: (Color online) Locations of nearest approach for both 
the Ef (dots) and E~ (solid) bands for the (a) x=0.12 and (b) 
£=0.19 cases with the color scale of the individual pseudogap 
weightings and . This illustrates the restructuring 
of the tight-binding Fermi surface in the antinodal direction, 
shown as the black dots of the curve in (b), related to the Eg 
band. 



Dirac points. For angles away from the node, there are 
no E~ = points. It is for these angles that the minima 
in energy serves as an effective Fermi surface. 

When these nearest approach energies are traced in k- 
space one obtains the inset of Fig. [1] demonstrating the 
Fermi pockets which are an essential part of the YRZ 
model. As has been previously described, the and 
E~ bands have relative weightings W£ and Wj^, given 
in Eq. ([4| which are dependent upon both the sign of the 
bare dispersion, as well as modified by the presence of the 
pseudogap, such that Wj[ + = 1. These weightings 
are incorporated into the nearest approach momentum 
contours of Fig. [2] through a color scale as shown by two 
vertical columns with black equal to one and orange equal 
to zero. 

Of course, each branch will have, possibly, separate 
nearest approach surfaces. Examples of these are shown 
in Fig. [21 In Fig. 2(a) which contains a large pseudogap, 
the E 



band, shown as a continuous curve, dominates in 
weighting along the frontside of the Fermi pocket closest 
to (0,0), forming a strong Fermi arc, shown in black (the 
top of our color scale). The backside of the Fermi pocket 
closest to the AFBZ boundary has a small weighting due 
to proximity to the pseudogap located there and is shown 
in orange which is the bottom of our color scale. In a 
recent ARPES experiment^ the backside of the pocket 
has finally been resolved and is indeed seen with much 
lower intensity as compared with the front side. The sec- 
ond momentum space contour, shown as solid dots, in 
Fig. 2(a) are associated with the branch. For most 
of this curve the intensity falls in the middle of the range 
shown, which is red on our color scale. Note that while 
these nearest approach contours are well defined, Fig. [2] 
tells us nothing about how close to zero the approach 
energies might be. This can be very important in de- 
termining how effective a process involving such energies 
might be. Later we will return to this issue. 

Fig. |2(b)| shows how the closest approach surfaces 
change as the pseudogap becomes small. Here for dop- 
ing x = 0.19, which is only slightly less than x c — 0.20, 



the Ef band has substantial weight near the antinodal 
'wings', shown as black dots. One can see the reformation 
of the tight-binding Fermi surface as the combination of 
Ef and E~ black regions. The side of the pocket nearest 
to the AFBZ boundary now has a more uniform weight- 
ing, approaching zero in value for x = x c . The contours 
also take on the shape of the AFBZ with which they al- 
most overlap. In the antinodal direction the Ej band 
has two nearest approach values as does E~ in the Fermi 
pocket region. 

We can also examine the size of the nearest ap- 
proach energies as a function of angle, along the strongly 
weighted region of Fig. 2(a) This is shown in Fig. O 



where we compare the cases with either a pseudogap or 
a superconducting gap to the case with both. The sur- 
face traced by minima in the E~ band coincides with 
the Luttinger pocket and actual E~ = Fermi sur- 
face in the absence of superconductivity, which is im- 
portant for correctly describing low energy nodal exci- 
tations. The strongly weighted part of the Ef band in 
the nodal direction occurs at a high energy relative to 
a zero energy Fermi surface that exists in the E~ band 
in the (7r,7r) direction. Comparison with ARPES data 
with energy scales relative to the Fermi level will there- 
fore be dominated by the E~ band. It is an important 
note that the nearest approach contours are not shifted 
when the superconductivity is switched on or off. We 
make two calculations of E~(0): one without the su- 
perconductivity in which we denote the energy of near- 
est approach by A na (T = T c ) (red short-dashed curve) 
and one with superconductivity denoted by A na (T = 0) 
(solid black curve). We can also evaluate directly the 
size of the superconducting gap on the nearest approach 
contour, A sc (9) and this is the long-dashed green curve. 
In the pseudogap case, for angles less than some crit- 



ical angle, yrz , 
pseudogap, is finite 



A na (T — T c ), which contains only a 



These angles 9 < 



represent 
In 



the region along the 'wings' of the Fermi surface, 
this sense 9 yrz marks the corner of the Fermi pocket. 
For 9 > Q yrz , A na (T — T c ) is uniformly zero along the 
strongly weighted side of the Fermi pocket. In the YRZ 
model, in spite of introducing a d-wave pseudogap over 
the entire Brillouin zone the net effect is that it only ap- 
pears dominant in the wings, and not in the pocket itself 
where the excitation spectrum shows no gap, as is also 
seen in ARPES data. Fig. [3] also includes A na (T = 0, 9) 
(solid black curve) which has both gaps present. The 
value of A na (T = 0,6) is dominated at 9 < 9 yrz by 
the pseudogap contribution, and by the superconduct- 
ing contribution for 9 > 9 yrz . However, there appears 
here a natural energy scale which is well approximated 
by A na (T = 0,6) « A tot (9) = JKJ C (9) + AjjT^ff), 
the root sum of squares of the superconducting and pseu- 
dogap contributions. This square root is represented by 
the open blue circles which overlap precisely with the 
A no (T = 0) curve. 

Comparison with experiment^ yields a complication 
in that there appears to be minimal difference between 
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A na (T c ,9) and A na (T = 0,8) for small angles. In order 
to obtain a similar feature that A na (0,8) « A na (T c , 0) 
for small angles, corresponding to the region off of the 
Fermi pocket, one would infer from the equation for 
Atot ($) that the new superconducting gap should be zero 
off of the Fermi pocket and be equal to A sc (9) from Fig. [3] 
for large angles, corresponding to the region on the Fermi 
pocket. Guided by Fig. [31 we have chosen a A sc (9) which 
is zero for small angles and smoothly transitions to be 
equivalent to its value in Fig. [3J for large angles. Fig. Q] 
demonstrates the result of such a constraint on A sc (9) 
(long-dashed green curve) as well as the expected over- 
all gap A na (0,9) = y/ A* c {0) + A* a (T c , 8) (solid black). 
This procedure results in a superconducting gap profile 
which appears either highly nonmonotonic, or that is only 
present on the Fermi pocket for 9 > 9 yrz , dropping off 
quickly beyond that. One can apply this simple analy- 
sis to the experimental data taken from RefJ^ which we 
have shown in the inset of Fig. |4j We see that the su- 
perconducting gap (solid green curve) is nonzero only in 
the nodal region up to approximately 23°, i.e. between 
23°-45° 



The issue of how the superconducting gap presents 
itself on the Fermi surface in complicated electronic 
systems is important. Even in conventional metals, 
for which the electron-phonon interaction is responsi- 
ble for the superconducting condensation, the gap is 
found to be highly anisotropic and does not exist for 
certain solid angles in momentum space where there is 
no Fermi surface because it is gapped out by the crys- 
tal potentiaL 23 i 24 ' 25 In the cuprates, significant devia- 
tions from the simplest d-wave gap, Eq. JT]), have been 
known for some time in Bi2Sr2CaCu20s+^^ and are 
seen in ARPES data which resolves bilayer splitting in 
(Pb,Bi) 2 Sn2CaCu08+,5>2£ A nonmonotonic gap was also 
observed in Nd2- a; Ce a: Cu04 ) 2 ^ an electron doped super- 
conductor. This is also expected in theoretical models 
with pairing based on spin nuctuationsj 29 i 30 ' 31 i 32 Numer- 
ical solutions of the gap equation in such cases show that 
many higher d-wave harmonics are present in the gap 
function. When a pseudogap is present, as is the case in 
the underdoped cuprates, one expects it to prevent the 
lowest order harmonic contribution to the superconduct- 
ing gap from having its full amplitude. This amplitude 
could even be forced to zero as is indicated in the data 
shown in the inset of Fig. U for the regions away from the 
Fermi pocket. The Fermi pocket is the only region where 
the normal state supports zero energy excitations. 
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FIG. 3: (Color online) energies of nearest approach, 
A na , versus angle, 8, for the x = 0.12 doping case: with 
superconductivity (solid black line) and without supercon- 
ductivity (short-dashed red line). The long-dashed green line 
is A sc (8) as described in the text. Marked with open cir- 
cles is the robust approximation that A na (T = 0,9) « Atot, 
where A tot (6) = \/AI c (6>) + Al a (T c , 8). The angle at which 
A na (T c ) = marks the edge of the Fermi pocket and thus 
defines 9 yrz , as measured from (ir,n) and described in the 
text. 



.ARC/ 



i n9 

A pg cos 



{ A9-n/2) 



<0 o ) 



> 



This Ap^- C (9) is essentially a shrunken cos(2(?) which 
goes to zero at 9q rather than 9 = 45°. For refer- 
ence, the open blue circles in Fig. 2] show Eq. (fT!?)) for 
9 = 9 yrz . This arc model pseudogap profile shows mini- 
mal difference from what would be extracted on the heav- 
ily weighted side of the pocket in the YRZ model. In this 
respect the two models are equivalent, the arc is just one 
side of the pocket. 



Of additional note, the angular profile of the pseudo- 
gap, A„ a (T c , 9) can be compared to another model, such 
as the arc model, which has previously been used to fit 
electronic specific heat dat a 11 ' 33 as well as understand the 
existence of two energy scales in Raman spectra.— The 
arc model, named as such because it produces a Fermi 
arc rather than a Fermi pocket, places a pseudogap over 
only a portion of the large tight-binding Fermi surface, 



IV. RAMAN IN THE YRZ MODEL 

Experimental Raman spectra suffer from a loss of the 
Bi s signal with increased underdoping as well as a strong 
signal from inelastic scattering at higher frequencies and 
complications due to surface effects. Still, given these 
issues, work has been done which is able to see clearly 



6 



0.25 



0.15 



0.05 




20 30 40 

(deg) 



50 



FIG. 4: (Color online) Nearest approach energies as a function 
of angle, 9 [measured from (n, n)], at T = (solid black line) 
and T = T c (dashed red line) for doping x — 0.12. Here the 
superconducting gap profile A sc (8) (long-dashed green line) 
is chosen so that A na (T = 0,9) « A na (T = T c ,0) for small 
9, which is different from Fig. [3] and is governed by A sc (8) — 
\f A^ a (0, 8) — A'n a (T c , 9). Inset shows a similar analysis of ex- 
perimental data from Kondo et ai^ for an underdoped sample 
measured at T = 10K (low temperature) and T — 40K (above 
T c ), where we take A sc (8) = y/A% a (WK,9) - Ai a (40K,9). 



two energy scales; 15 i 35 i 36 one in the antinodal B lg sig- 
nal which increases and the other in the nodal B2 ff signal 
which decreases with underdoping. While ARPES exper- 
iments are able to resolve fc-space dependence, Raman 
spectra sample a finite region of fc-space associated with 
nodal or antinodal polarizations. These regions may even 
overlap somewhat. We seek to illustrate how some of the 
intricate details seen in ARPES, namely the Fermi sur- 
face restructuring and the superconducting gap profile, 
impact upon Raman spectra. 

Fig. [5] shows a sample calculation for the £=0.10 case. 
Each frame contains Raman spectra for the pseudo- 
gap case only (dashed green), superconducting case only 
(dash-dotted red) and both (solid black). In the top 
frame the Bi g peak, 2 A an, is marked. For this dop- 
ing, the 2A J 4 J y peak occurs at a higher frequency than 
both the 2A pg and 2A SC points due to the additive im- 
pact of these two gaps on the antinodal gap. In the A sc 
only curve, the peak occurs just before the 2A SC point. 
In the simplest BCS calculations based on free electron 
bands, the position of this peak would coincide exactly 
with the value of twice the superconducting gap ampli- 
tude. However, in more complicated band structure mod- 
els, as we are using here, and with a superconducting gap 



defined over the entire brillouin zone rather than just on 
the Fermi surface this need not be the case. In addi- 
tion, the Raman cross-section can display other features 
such as the second peak seen in the red dash-dotted curve 
of the top frame of Fig. [5] which have their origin in the 
band structure. In the bottom frame, the B2 9 peak, 2Ajy 
(first peak in the B2 S spectrum), is marked on the solid 
black line, and occurs at lower energy than both 2A SC 
and 2A pg . There are two main factors which contribute 
to this lower value. First, the value of a gap in the B2 S 
spectra in standard BCS theory falls not at its peak but 
beyond the peak, at a point of inflection. Additionally, 
since the B 2s vertex samples largely the nodal direction, 
it may not contain the maximum input gap values which 
occur in the antinodal direction, depending upon the ac- 
tual size and shape of the Fermi pocket. 

To elaborate, as mentioned in Fig. [3j the edge of the 
Fermi pocket is marked by an angle 9 yrz which changes 
quite drastically with doping. If we instead measure this 
angle from (0,0), which we will distinguish by calling 
(shown in the upper right inset of Fig. [6]), we 



it 



^yrz 



would find that 4> yrz has a large value approaching 45° 
at x = 0, and for dopings at or above the QCP, x c , 
has some small value defined by the angle from (0,0) to 
the point where the tight binding Fermi surface inter- 



sects the Brillouin zone boundary. As shown in Fig. 2(a) 



the edge of this pocket can be quite heavily weighted, 
and as a result any experimental probes that are con- 
cerned with excitations along the Fermi surface will see 
a maximum superconducting gap value at the edge of 
the Fermi pocket near the angle 4> yrz where a maxi- 
mum value coincides with a strong weighting. In a the- 
oretical calculation, the input of a superconducting gap 
with a maximum value Ag C (a;) [see Eq. (JTJ] such that 
A sc (a;, 4>) — Aj? c (a;) cos(2</>) will show features in the Ra- 
man spectra at a value A SCeff (x) = A sc (x, 4> yrz ). Fig. [6] 
demonstrates the size of this effect in the YRZ model. 
Shown in solid black is the input superconducting dome 
given by Eq. (J5J) while the open purple squares show A SCeff 
in the YRZ model. This is obtained by finding the an- 
gles <pyr Z (x) across the doping phase diagram. Although 
this effect is small in the absence of pseudogap; in the 
presence of pseudogap, this alone results in a factor of 2 
difference between the apparent gaps, as compared with 
the actual input gap, across the phase diagram. This is 
seen clearly in the upper left inset of Fig. [6j which shows 
just the value of cos(2(f> yrz ) versus doping. 

The shrinking of the Fermi pocket with underdoping 
has additional important effects on the Bi fl antinodal re- 
sponse. Fig. [7] shows an overlay of the effective Fermi 
surface onto the Bi ff vertex strength, ("f^ 19 ) 2 - The Bi s 
vertex cannot properly sample the Fermi surface near 
the edge of the Fermi pocket for strongly underdoped 
cases. Although this is an extreme example, illustrat- 
ing the large range of pocket size, this effect can be seen 
clearly over a much smaller doping range. Fig. [8] shows 
the YRZ Raman B l9 spectra with the doping-dependent 
Gutzwiller factors removed to better illustrate this ver- 
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co/t 

FIG. 5: (Color online) Raman response curves for the a;=0.10 
case. The top and bottom frames show the Bi 9 and B2 9 
responses, respectively. The antinodal gap, 2Aajv, and the 
nodal gap, 2Ajv, are defined as the dominant peaks in their 
respective Bi 9 and B2 S curves. As well, the actual input gap 
values 2A SC = 0.2to and 2A P9 = 0.6to are marked. 



tex strength effect. For doping above x = 0.17, there 
are extremely strong low frequency peaks which mark 
the antinodal gap. For dopings below x — 0.17 the low 
frequency signal has become extremely suppressed. The 
horizontal dashed line helps the reader to see the two dis- 
tinct regions of doping with the boundary between the 
regions occuring at x = 0.17. This doping falls consider- 
ably below its critical value, £ c =0.20, which defines the 
quantum critical point at which the pseudogap sets in. 
In both regimes, the Raman vertex favors the antinodal 
direction. While it can vary somewhat, both in shape 
and in amplitude, since the hopping parameters in the 
electronic dispersion curves depend on x, this is not the 
important effect that we wish to emphasize here. What 
is more important is the evolution in the topology of the 
Fermi surface, or more correctly the surface of closest ap- 
proach, that occurs for x < x c . As the pockets begin to 
shrink they move away from the antinodal point and the 
intensity of the Raman signal consequently drops due to 
an effectively smaller value of the vertex projected onto 
the pocket. To get a large reduction it is necessary that 
the distance in momentum space between the end of the 
pocket and the Brillouin zone boundary at the antinodc 
is not too small. Our calculations show that this occurs 
rather abruptly at a doping of x — 0.17. 

So far we have not accounted for damping. The in- 
clastic scattering in the cuprates, which is known to be 
large, can broaden the Raman response considerably and 
provide a background that remains in the Fermi liquid 
state d 7 ' 38 In our formalism it would correspond to the 



incoherent part of the Green's function which should add 
a second piece in Eq. (J3]) which involves only the coher- 
ent part. It has been shown that the Raman, optical, 
and quasiparticle scattering rates are all strongly depen- 
dent on frequenc y 38 ' 39 ! 40 ' 41 In the superconducting state, 
these scattering rates are reduced at low energy by the 
opening of a superconducting gap and there is therefore 
no justification for broadening the B2 9 spectrum which 
samples mainly the low energy nodal region. However, 
the antinodal Bi s spectrum contains peaks that generally 
fall at higher energies, particularly in the underdoped re- 
gion of the phase diagram, and we therefore expect rather 
large scattering rates and considerable broadening. For 
experimental purposes, one would need to know the pseu- 
dogapped normal-state Raman response to compare to 
the superconducting Raman spectra in order to see indi- 
cations of the strength of the superconducting state in the 
combined spectra. Fig. [5] gives such a comparison for the 
Bi g (right frame) and B2 S (left frame) spectra where the 
Bi g curves have been broadened with a doping dependent 
elastic scattering of T — 0.2t(x) while the B2 9 have been 
kept near the clean limit. There is an apparent shift of 
Big spectral weight, shown by the shaded yellow region, 
upon the addition of a superconducting energy scale. The 
differences between the normal and superconducting Big 
spectra become much less pronounced for increased un- 
derdoping. This behaviour occurs at all dopings where 
the pseudogap energy scale dominates, and the super- 
conductivity provides only a very small additional con- 
tribution to the antinodal Raman. We note also that 
the overall strength of the Raman signal decreases with 
decreasing doping as noted in experiments! 15 ' 36 As we 
have described before, this can be traced to a smaller 
overlap between the appropriate Raman vertex weight- 
ing factor in the Brillouin zone and the Fermi surface 
pocket. The case of the B2 9 response yields a different 
result. The region of interest occurs at a lower frequency, 
on the scale of the superconducting energy gap, while the 
effect of the pseudogap is most apparent only at higher 
frequency. However, the pseudogap does suppress the 
lower frequency normal-state Raman background. This 
makes the effect of the superconducting gap very appar- 
ent relative to this background even at small dopings. 
Note that in the lower left frame, the shaded yellow re- 
gion at the lowest doping remains substantial compared 
with its value in the upper frame. This is distinct from 
the case in the right hand frame where the shaded yellow 
region at low dopings has essentially disappeared. 

We have calculated the temperature evolution of these 
peaks for the case where the superconducting and pseu- 
dogaps are full <i-wave gaps over the Brillouin zone, given 
by Eqs. ([T]) and ^ and where we have assumed A° g to be 
constant with temperature and A° c to have a BCS tem- 
perature dependence. The top frame of Fig. [TU] shows 
that the Bi g peaks, normalized by their zero tempera- 
ture values, do not follow the input BCS temperature 
dependence for the superconducting gap, but follow in- 
stead a temperature dependence related to the combined 
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0.15 0.2 
Doping x 

FIG. 6: (Color online) Plot of the superconducting dome 
A sc (x) (black) and the reduced effective gap A SCeff (:r) = 
A sc (a;) cos(2^> Hrz ) (open purple squares). This effective gap 
represents the maximum gap value observed on the Fermi 
pocket, which occurs at an angle <j> yr z as measured from (0,0) 
(see right inset). The factor cos(20 yrz ) is plotted as a function 
of doping, x, in the left inset. 



magnitudes of A SC (T) and A pg . For the extremely un- 
derdoped cases, there is almost no discernible tempera- 
ture dependence. For increased doping, the pseudogap is 
reduced, and as x c is approached, the assumed BCS tem- 
perature dependence returns. This is distinct from the 
temperature dependence in the nodal direction which is 
dominated by the temperature dependence of the super- 
conducting gap as is shown in the bottom frame of Fig. 1101 
where it is compared to the input BCS A(T)/A(0). 

We also wish to consider a cutoff superconducting gap, 
a gap only present along the Fermi pocket, as suggested 
by ARPES experiments. As a first attempt, we can cut 
off the superconducting gap sharply at the edge of the 
Fermi pocket. Although this sharp cutoff would result in 
clear features in nearest approach energies, the sharpness 
should not create additional features in the Raman spec- 
tra due to the averaging effects of the Raman vertices. 

To state explicitly, when we state 'A sc cutoff', we are 
referring to 



-(cos k x a 



cos k y a) 



< 



< - 

^ 2 



0. 



otherwise. 



(20) 

With this A sc cutoff, the gaps in the nodal direction are 
largely unaffected, and the temperature dependence of 
the corresponding B2 9 peaks follows exactly as in Fig. [TOl 
(lower frame). This is not the case for the Bi g antinodal 
results, where the superconductivity has been removed. 
These results are shown in Fig. [TTJ The open blue tri- 
angles are for the x—0.16 case from Fig. [TUl which had 
no A sc cutoff. Upon inclusion of the cutoff, the results 
change drastically to the solid blue circles, which show no 
temperature dependence (due to the lack of temperature 




I 




(b) 








V 


^\ 











K 



(c) 



4.5 
4 

3.5 
3 

2.5 
2 

1.5 
1 

0.5 
() 



4.5 
4 

3.5 
3 

2.5 
2 

1.5 
1 

0.5 
() 



FIG. 7: (Color online) An overlay of the Fermi pocket and 
Fermi surface wings (in white) for the band, on a back- 
drop of the Bi s vertex strength, (7 Bl9 ) 2 , as in Eq. (fT5}, for 
dopings: (a) z=0.05 (b) a;=0.16 and (c) x=0.19. This il- 
lustrates an expected loss of Bi s signal in the underdoped 
region as the pocket, which dominates the superconducting 
state, shrinks away from regions of strong (7 s19 ) 2 , in addi- 
tion to loss associated with the considerable weakening of the 
vertex itself. 



dependence assumed for the pseudogap). Although not 
shown, all dopings below x=0.16 show the exact same 
result that the dominant peak is purely pseudogap and 
independent of T . However, as the doping is increased 
above a;=0.16, there is reconstruction of the Fermi pocket 
in the antinodal wings, shown in the inset of Fig. [IT] for 
the a:=0.18 case. For these cases, the superconducting 
gap is present in the antinodal direction, and as a result, 
the temperature dependence of the Bi s peaks returns. 
This change in Fermi surface in the antinodal direction 
occurs over a relatively small doping change, Aa;=0.01- 




B2g Nodal Region 



FIG. 8: (Color online) Raman Bi 9 spectra versus frequency 
for a range of doping values wherein the gl (x) prefactor has 
been removed. Below x—0.17, the Xs ls amplitude is small at 
low frequencies but maintains a roughly constant peak ampli- 
tude. Above £=0.17, we see a strong increase in low frequency 
amplitude caused by the reformation of Fermi surface in the 
antinodal region which is of interest for the Bi 3 vertex. A 
dashed horizontal line has been added to help the reader dis- 
tinguish these two regions. 



0.02, and results in drastic modification to the temper- 
ature dependence of the B i£( Raman spectra. This gap 
cutoff could explain the recent data of Guyard et al.^ 
which sees a large change in the Bi g temperature depen- 
dence over a small range of doping. Furthermore, dop- 
ings which present little or no temperature dependence 



in the B ig spectra, despite a strong B2 9 nodal gap A^v, 
are evidence of a pseudogap suppressing the antinodal 
part of the Fermi surface. Experimental observation of 
this lack of temperature dependence continuing below T c 
is evidence to support the idea that the pseudogap line 
continues through to T=0 in the phase diagram, and does 
not terminate at the top of the dome. 

Comparing Raman spectra across the phase diagram 
is complicated due to doping dependent changes to the 
vertex weightings, as well as modifications to the Fermi 
surface. To further emphasize this point, we calcu- 
late the low frequency slope of the nodal Raman re- 

sponse, S B 2 9 = 



, across the doping phase di- 
agram. In BCS theory of d-wave superconductors, this 
low frequency dependence goes linear with frequency as 
XB2g ^ W /Asc- Since for most dopings there is little 
impact due to the pseudogap on the nodal direction re- 
sponse, the effects of the vertex component, as well as 
Fermi pocket size, can be seen clearly. Fig. [12] shows 
these slopes versus doping as open red circles with the 
solid red line a guide to the eye. Taken directly, there 
is a factor of 5 difference across the region shown. Wc 
can largely remove the doping dependent vertex factors 
by scaling out [t 1 (x)g t (x)] 2 , and can compensate for the 
superconducting gap by multiplying by A sc (x). This re- 
sults in the dashed blue curve in the upper right inset of 
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FIG. 9: (Colour online) Raman response in the YRZ model 
for the B2 9 (left frames) and Bi g (right frames) polariza- 
tions. The dashed red curves are for the normal state, and 
act as a background for reference. In the underdoped region 
the pseudogap background dominates, and superconductivity 
plays little role in the Bi ff spectra. This is emphasized by the 
shaded region between the curves above the background peak. 
The effect of A sc on the B2 9 curves is clear, even at low dop- 
ings, by a strong B2 9 signal relative to the background. Also 
note that the y-axis scales are modified with doping, which is 
an effect of the gi(x) prefactors. 



Fig. [T2l With the factors removed there remains a dop- 
ing dependence, which varies by a factor of 2 across the 
doping range shown. To compensate for the effect of the 
pocket size, we replace the A sc (a;) factor with a A SC(stf (aj) 
factor, which is shown in the top left inset of Fig. [121 To 
within 5% in the values here, the simple BCS low fre- 
quency scaling is maintained since this is purely nodal 
physics and should be relatively unaffected by the higher 
energy pseudogap scale. This does illustrate however, 
that the indirect effect of the pseudogap of modifying 
the Fermi surface has an important and strong impact 
on the low frequency B2 9 spectra. This shows that it is 
the value of the superconducting gap on the Fermi pocket 
that is important, and not its value at other points in the 
Brillouin zone. 
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FIG. 10: (Color online) The normalized A A n/A A n{T = 0) 
(top frame) and Ajv/Ajv(T = 0) (bottom frame), as would 
be extracted from spectra similar to those in Fig. [5] but taken 
at various finite temperature, are plotted versus T/T c for a 
range of dopings. The nodal gap, Ajv, shows BCS-like tem- 
perature dependence across all dopings, while the Aan di- 
verges from BCS temperature dependence in the presence 
of a modest pseudogap. Although not shown, this diver- 
gence from BCS behaviour is reasonably well-described by 
A(T) = ^AUT) + A* g . 



Finally, we wish to examine the progression of zero 
temperature peaks in the Raman response as the doping 
is varied over the phase diagram. Fig. [13] shows such a 
plot, which includes the input 2A sc (x) dome (solid black 
curve) and 2A pg (a:) line (dash-dotted green line). Also 
included is 2A SCcff (x) (open purple squares), the maxi- 
mum gap value on the Fermi pocket, from Fig. [5] All 
zero temperature calculations are for a superconducting 
gap which exists over the entire Fermi surface and not 
just on the pocket, as in Eq. (JTJ). The open red dia- 
monds mark Bi g peaks, while open black circles mark 
B2 S peaks. One first notes that the B2 9 peaks lie sub- 
stantially below the input 2A SC dome and follow a curved 
shape much like the 2A SCoff (x) modified dome. In prac- 
tice, the 2A SCoff point should always lie beyond the B2 9 
peak at a point of inflection (as mentioned in discussion 
of Fig. [SJ). Here, if one were to divide 2A SCcff (x) by the 
B2 9 peak values, you would find a roughly constant value 
of w 1.4. The Bi s peaks show additional features. For 
dopings above £ c =0.20, the B ls peaks fall exactly on the 
2A SCoff dome. Just below x c , for dopings x=0.17— ►0.19, 



FIG. 11: (Color online) Normalized antinodal peak frequen- 
cies as a function of temperature. Here, solid points illustrate 
the effect of cutting off superconductivity in the antinodal 
regions as in Eq. (|20p . This cutoff would result in no tem- 
perature dependence in the Bi 9 peaks for any doping below 
£=0.16. Above x=0.16, the Fermi surface begins to restruc- 
ture in the antinodal directions and begins to be sampled by 
the Bi 9 vertex. Once this occurs, BCS-like temperature de- 
pendence returns. The inset shows the reformation of the 
antinodal Fermi surface for the £=0.18 case. 



the Bi s peaks lie just above the 2A SCoff dome, due to the 
presence of a small pseudogap. For these dopings, there 
is still Fermi surface in the antinodal region in the form 
of small 'loops' which are nearest approach surfaces for 
the band, as is illustrated in the inset of Fig. [TTJ for 
the £=0.18 case. There is also a jump in this curve be- 
tween £=0.16 and £=0.17. For dopings above £=0.17, 
the Big spectral peaks coincide with the E+ energy on 
the nearest approach surface in the antinodal direction. 
For dopings below £=0.17, the Ef 'loops' disappear al- 
together. This results in a shift from the E+ band to the 
E~ band, which happens to have a higher energy of near- 
est approach in the antinodal direction. As a result, the 
Big peaks jump up in energy as x is changed from 0.17 to 
0.16. We will refer to this jump at £=0.16 as being at the 
onset doping, x onse t , the doping below which the system 
is dominated by pseudogap-created Fermi pockets. 

We now have the image of three distinct regions: 
% < x onset , which is dominated by the pseudogap and 
the Fermi pockets; x onS et < x < x c , which is a region 
wherein strong Fermi surface reconstruction occurs; and 
finally x > x c , where the pseudogap is zero and the sys- 
tem is dominated by the A sc energy scale with a large 
tight-binding Fermi surface described by a Fermi liquid 
model. Of these three regions, the one that is most com- 
plicated to understand quantitatively is the middle range, 
where the topology of the Fermi surface becomes particu- 
larly complex and includes Fermi pockets as well as other- 
pieces of Fermi surface with complex geometry around 
the antinodal points as seen in the inset of Fig. [TT] It 
is however precisely this region where the Big Raman 
scattering cross-section shows the most rapid variation, 
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FIG. 12: (Color online) Low frequency slope of the B2 9 Ra- 
man response, SB2 g = ^— ^p-J , versus doping. In the 

top right inset, doping dependent factors [t' (x)gt(x)] 2 have 
been removed to attempt to verify the low frequency scaling 
°f XB2g K w/A ac . The top left inset shows that the scaling 
relation is changed to Xs2 9 tx cj/A SCoff in order to compare 
slopes across all dopings. To guide the eye in the top left in- 
set, we have included a horizontal line at the average value of 
the points and included error bars of ±5% of the point values. 
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FIG. 13: (Color online) Gaps versus doping. Shown are in- 
put values of 2A sc (a:) (solid black) and 2A pg (a;) (dash-dotted 
green) given by Eqs. ([5jl and ©, respectively. Open red dia- 
monds are the Bi 9 peaks extracted at zero temperature. Open 
black circles are the B2 9 peaks at zero temperature. The 
open purple squares give the effective superconducting gap, 
2A S c eff (x), as described in Fig. [6] In the overdoped region, 
the Bi 9 peaks fall on the A SCrff (x) line, while the B2 9 peaks 
fall well below. 



as seen in Fig. [5J and has the greatest promise for future 
experimental investigation. A second interesting feature 
to note about the phase diagram traced in Fig. [13] is 
that the values of the position of the B lg Raman peaks 
(open red diamonds) fall below the input pseudogap val- 
ues (green dash-dotted curve) in the highly underdoped 
region near the lower end of the superconducting dome. 
These reduced values coincide with the value of the E~ 
band on the antinodal nearest approach surface. The 
physics behind this reduction in peak position relates 
to the shape of the Bi g Raman vertex and its overlap 
with the Fermi pocket as is seen in Fig. 7(a) The Fermi 
pockets have shrunken so much that the peak in the Bi s 
Raman cross-section comes from sampling a momentum 
region which is displaced from the antinodal direction 
towards the nodal direction and in this region, the pseu- 
dogap is reduced from its full amplitude A° (x), by the 
modulating factor [cos(fe x a) — cos(k y a)) of Eq. ([2]). 

Fig. Q3] contains a comparison of the results of our cal- 
culations to the experimental data of Guyard et at, taken 
from Ref^. The top frame contains two y-axis scales. 
The left scale is for our calculation, with dimensionless 
quantity A/to an d applies to the open points, while the 
right hand scale has units of meV and applies to the 
experimental data marked with filled points. The ex- 
perimental and theoretical scales are shown with equiv- 
alent values for the assumed to = 175 meV in the YRZ 
model. It is clear on first glance that the energy scale for 
our calculations agrees well with experiment despite not 
being rigorously fit to this specific experimental system 
of Hg-1201. The experimental data sets show that Ajv 
(solid black upward pointing triangles) and Aan (solid 



red downward pointing triangles) have very similar val- 
ues above some doping around x « 0.16. This feature is 
not captured in our calculation for the An (open black 
circles) and Aan (open red diamonds), as Ajy main- 
tains a lower value than Aan for all dopings. We note 
however that A SCe/ , (x) (open purple squares) has simi- 
lar values to Aan (open red diamonds) in the overdoped 
region. This illustrates that the actual An doping depen- 
dence, as seen in experiment, may not necessarily follow 
a straightforward superconducting dome, but rather an 
effective dome, given by A SCoff (x), which includes Fermi 
surface restructuring effects due to the pseudogap. Our 
parameters could be improved by scaling down our en- 
ergies, which could equivalently correspond to to being 
smaller by a factor of w 1.1. 

In the lower frame of Fig. 1 141 we plot the experimental 
ratio Aan/ An as solid blue triangles. We can clearly 
see two doping dependent regions of behaviour in this 
data: x < x onse t and x > x onse t, where x onS et marks the 
change in slope around £=0.16. We seek to distinguish 
the difference between x onse t and x c in this data. The 
antinodal to nodal gap ratio, Aan / An (open red circles) 
as well as Aan / A SCe£f (shown in open black squares) are 
shown for the YRZ model. The Aan/ An calculations 
also show two distinct regions: x < x onse t and x > x onS et, 
where x onse t=Q-17 ■ For dopings above x onS et, the Aan 
and An curves do not have the same value, and instead, 
have a roughly constant ratio of w 1.4. When we instead 
normalize Aan by A SCcff , resulting in a constant ratio 
of 1 for x > 0.16, we find better agreement with the 
experimental observation of the Aan/ An ratio. 

An important point is that all three data sets (2 the- 
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FIG. 14: (Color online) Top frame: The data of Guyard et 
al. [Ref.— Fig. 2(a)] for Aajv and Ajv is shown in solid 
points (right hand axis) while the YRZ results for Aajv, Ajv 
and A SCoff are shown in open points (left hand axis). Bot- 
tom frame: Comparison of the antinodal to nodal gap ra- 
tio Aajv/ Ajv for experiment (solid blue triangles) and YRZ 
model (open red circles) as well as to the YRZ result for 
Aajv / A S c off (open black squares) . 



oretical, 1 experimental) in the lower frame of Fig. 1141 
show only two distinct regions of doping dependent be- 
haviour: x < x nset and x > x nset- There is virtually 
no signature of x c in this data. We can analyze this fur- 
ther, since the x onS et marks the dominance of A pg over 
A SCoff in the phase diagram of Fig. [T3J all x < x onse t 
are dominated by the pseudogap. If we assume that the 
pscudogap energy scale maintains linearity, we can trace 
a simple linear fit of the points just below x onse t and 
extend that line above x onS et ■ This will result in an esti- 
mate of the location of the doping at which the pseudogap 
vanishes, and hence, the doping where one might look for 
evidence of a zero temperature QCP. Indeed this extrap- 
olation gives x c « 0.2, as was input into the theory, and 
it provides an estimate from experiment that x c ks 0.19. 
It is with this simple interpretation that we reestablish 
the existence of three regions of importance in this data: 
x < x onset (where A pg > A SCotf ), a; „ set < x < x c (where 
Apg < A SCoff ) and x > x c (where A pg = 0). 



V. CONCLUSIONS 

The YRZ model provides a formalism whereby a pseu- 
dogap opening up about the AFBZ boundary recon- 
structs a large Fermi surface into small Fermi pockets as 
the doping progresses towards the Mott insulating state. 
Comparison of the YRZ results to recent Raman spectra 
on Hg-1201 gives excellent qualitative agreement, and es- 
tablishes that its quantum critical point, associated with 
pscudogap formation, falls at doping £=0.19 inside the 
superconducting dome and is consistent with a zero tem- 
perature transition. We have shown that if the supercon- 
ducting gap is only present in the Fermi pocket region 
(nodal direction) then the B ls peaks show no tempera- 
ture dependence. However, beyond some critical doping, 
the Fermi surface reconstructs in the antinodal region 
of the Brillouin zone resulting in a large change in the 
Bi g peak temperature dependence over a small change 
in doping. We also understand the loss of B lg spectral 
amplitude in the superconducting state to be an effect 
of the presence of a dominant normal state background 
in the antinodal direction which grows in strength with 
underdoping. 

Tracing the major peak in each of the Bi g and B2 9 
response curves as a function of doping creates an appar- 
ent phase diagram in the YRZ model which illustrates 
the indirect impact of the presence of the Fermi pocket. 
The effect is that the Raman spectra can only sample the 
largest value of the gap on the Fermi pocket which is at 
the edge, towards the antinodal direction. This results in 
an effective superconducting dome A SCef( (ai), which can 
be significantly smaller than the input gap. This occurs 
for dopings where the pseudogap is strong, resulting in 
Fermi pockets which become smaller on approach to the 
antiferromagnetic state. We have also analyzed the low 
frequency slope of the B2 9 polarization, taking into ac- 
count the Fermi surface restructuring (pocket size) effect 
by replacing A sc with the effective gap, A SCoff . 

Through this work we have identified three regions 
of doping dependence in the phase diagram: the un- 
derdoped pseudogap-dominated region, the overdoped 
superconductivity-dominated Feremi liquid region, and 
the intermediate mixed region wherein the presence of 
pseudogap begins to erode the Fermi surface. In the 
case of Raman spectra, this intermediate region may be 
seen best in the temperature dependence of the antinodal 
peaks, but also as a sudden jump of the antinodal peak 
frequency with doping. 

Although our results are largely concerned with the 
formation of Fermi pockets, we understand that there 
are similar models of the pseudogap state which result in 
the formation of Fermi arc segments rather than Fermi 
pockets. The qualitative results given here will remain 
unchanged, as they require: 1) some disappearance of 
Fermi surface in the antinodal direction due to the pres- 
ence of a pseudogap, 2) a dominant superconducting gap 
value given by the maximum value on the Fermi surface. 
These two requirements should still be present in models 
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which result in arcs rather than pockets. 
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